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The results of difluse neutron scattering experiment on PbMgi/3Nb2/303 single crystal above the 
Burns temperature are reported. It is shown that the high temperature elastic diffuse component is 
highly anisotropic in low-symmetry Brillouin zones and this anisotropy can be described using Huang 
scattering formalism assuming that the scattering originates from mesoscopic lattice deformations 
due to elastic defects. The qualitative agreement between this model and the experimental data 
is achieved with simple isotropic defects. It is demonstrated that weak satellite maxima near the 
Bragg reflections can be interpreted as the finite resolution effect. 

PACS numbers: 61.72.Dd, 77.84.Dy 



I. INTRODUCTION 

Ferroelectrics with diffuse phase transition were first 
discovered by Smolensky et alfi more than fifty years ago. 
Now they are usually referenced as relaxor ferroelectrics^.. 
These materials demonstrate exceptional piezoelectric 
capabilities^ and very high dielectric constants in a broad 
temperature range making them promising materials for 
applications in electronic and electromechanical devices. 
The most distinguished features of relaxors are broad 
frequency dependent maximum of dielectric permittivity 
e(T)^ and the nonergodicity of low-temperature phased. 
Above the characteristic Burns temperature relax- 
ors behave similarly to normal displacive ferroelectrics 
- e{T) follows Curie- Weiss lawl and the soft mode is 
observed^.. Typical relaxors are perovskite-like mixed 
crystals with one or two sublattices randomly occupied 
by nonisovalent ions. The most studied relaxor com- 
pound is lead magnoniobate PbMgi/3Nb2/303 (PMN) 
having B-sublattice populated by Mg^+ and Nb^+ in 
proportion 1/2. Diffraction studies show that without 
application of strong electric field PMN crystals retain 
cubic structure at least down to to 5 K and no ferro- 
electric phase transition takes place^'^°. However it is 
clear that at least below Td « 650 K local structural dis- 
tortions develop ^ These distortions are refiected by 
quasi-elastic diffuse scattering (DS) that appears below 
Td in the vicinity of the Bragg refiections^'^'^. The scat- 
tering is highly anisotropic in reciprocal space and has 
different shapes in different Brillouin zonea^^. On cool- 



ing down below Td the intensity of DS monotonically 
increases and at cryogenic temperatures becomes com- 
parable with the intensity of Bragg reflections. There 
were proposed numerous interpretations of DS in relax- 
ors including scattering on highly anisotropic optic^"* and 
acoustioi^ phonon modes, scattering on static correlated 
ionic displacements in polar nanoregionslS and scatter- 
ing on static mesoscopic lattice deformations caused by 
polar defects^-. The Monte-Carloi^ii^ and molecualr 
dynamics^*^ simulations were also productively emoloyed 
in the recent years to describe the spatial distribution of 
the DS in relaxors. The temperature-dependent DS can 
be considered as one of the relaxor footprints but despite 
extensive studies its microscopic origin is not yet clear. 

Along with temperature-dependent DS a nearly 
temperature-independent components were observed in 
PMN and related compounds^^"^"^. First observation was 
reported in Refl2ll where the DS was studied in PMN 
using high-resolution synchrotron X-ray scattering. It 
was emphasized that the longitudinal (or radial) scans 
contain a component that does not change at least in 
80 K < T < 270 K range (Fig. [T]). This component 
was interpreted by the authors as Huang scattering due 
to elastic lattice deformations. In 2004 using neutron 
scattering Hiraka et al^ clearly observed weak temper- 
ature independent component of DS in (100) and (110) 
BZs both above and below Td and attributed it to the 
correlated displacements of Pb^+ ions due to short-range 
chemical order which is believed to be temperature in- 
sensitive up to T « 1000 K. This interpretation was sup- 
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FIG. 1: (Color online) The longitudinal (radial) scans of syn- 
chrotron X-ray diffuse scattering in PMN in the vicinity of 
(300) Bragg reflection for r=120 K (squares) and r=270 K 
(circles)^. 



ported in further study by Gehring et but opposed 
in molecular dynamics study of relaxor PSN^S where the 
radial DS was attributed to the chemical disorder. To 
clarify the nature of the temperature independent radial 
component of DS in lead magnoniobate we extended the 
measurements of neutron diffuse scattering to BZs of low 
symmetry. In this contribution we show that our data 
and the data from Reft2^ can be adequately described 
in terms of Huang scattering due to inhomogenous lat- 
tice deformations which are definitely expected in mixed 
crystals. 



II. EXPERIMENTAL DETAILS 

PMN single crystal was supplied by the Institute of 
Physics, Southern Federal University (Rostov-on-Don, 
Russia). The measurements were carried out at 5G- 
PONTA three-axis spectrometer installed in the reac- 
tor hall of the JRR-3M reactor (Tokai, Japan). The 
sample was wrapped in the Al foil and mounted in 
the furnace with (001) axis vertical. The experimental 
setup allowed heating up to 650K . Two-dimensional dis- 
tributions of the diffuse scattering intensity were mea- 
sured at T = 650 K in (300), (310) and (210) Bril- 
louin zones in quasielastic mode with fixed wavevector 
Ki = Kf ~ 3.84A^^ (30.5 meV). Energy resolution was 
of the order of 1 meV and Q-resolution (FWHM) was 
AQl = O.O6A-1, AQt = O.O2A-1, /^Qv = O.4A-1 
for the component along the scattering vector, in-plane 
transverse component and out-of-plane (vertical) trans- 



verse component respectively. Due to strong contam- 
ination by Al powder scattettering in (300) and (310) 
Brillouin zones some parts of 2-d maps on Fig. [5] are 
omitted. In order to clearly distinguish between Al pow- 
der scattering and PMN diffuse scattering the positions 
of Al powder coils are marked by dashed lines. 

III. RESULTS AND MODELLING 

In general the results of our experiment are consistent 
with the previous ones^^: in high-symmetry (300) Bril- 
louin zone the scattering is oriented along the direction 
of reciprocal lattice vector T which is denoted by the 
arrows on Fig. [2l and the longitudinal profile contains 
weak maximum on approximately 0.15 r.l.u. from the 
zone center. In zones of lower symmetry the situation is 
somewhat different - the dumbbell-like diffuse scattering 
shapes are not purely longitudinal. The areas with max- 
imal intensity are clearly shifted from the line Lr that 
passes through the reciprocal lattice point and is oriented 
along T . Even more clearly it is seen on 1-d slices of 2-d 
maps presented on Fig. [31 The slices are made along the 
straight lines in reciprocal space determined by the vec- 
tor equation ((^ — Co) •''') = where ^ = (/i, k, 0) and the 
center of the slice lies on Lt . These lines are shown in 
sohd on Fig. [51 

We look for the origin of the observed scattering in 
elastic lattice deformations that are unavoidable in mixed 
crystals. This kind of scattering was heavily studied in 
alloys and other systems with imperfect lattice. In the 
book of M. A. Krivoglaa2^ a useful formalizm for calcu- 
lating of the distribution of static ionic displacements in 
the lattice due to elastic defects is presented. The char- 
acteristics of the spatial distribution of atomic displace- 
ments due to elastic deformations are determined by the 
types of defects and by elastic constants of the crystal. 
In our analysis we used the most straightforward way of 
applying this formalism by assuming simple cubic sym- 
metry defects in PMN lattice whose elastic constants are 
known in advance^^. The Fourier components of lattice 
deformation field are described by the expression 

^m(q) = -NdC-Mci)pa.jini (1) 
q 

where Nd is the number of defects, q - reduced wavevec- 
tor, n — c{/q. The tensors Cy (q) = cajmniUm and 
Paij = CijimLaim are Calculated on the basis of elastic 
constants of the crystal Cijim- The elastic deformation 
tensor Laji describes the symmetry and " force" of defect 
with orientation a. The intensity of diffuse scattering 
due to the elastic deformations ([T]) is determined by ex- 
pression 

/H(Q) = |F(r)p^(Q.A„(q))2 (2) 

Q=l 

The scattering vector Q is a sum of reciprocal lattice 
vector T and reduced vawevector q, F[t) is the elastic 
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FIG. 2: (Color online) Isointensity contours of difTuse scattering in PMN at r=650 K in three BZs (top) and corresponding 
isointensity contours produced by model calculations described in text (bottom). The arrows through the zone centers denote 
the direction of reciprocal lattice vector T. The positions of Al powder coils are marked by dashed lines. The straight solid 
lines perpendicula to vector T mark the positions in which the slices depicted on Fig. |3]are made. 
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FIG. 3: The 1-d slices of 2-d diffuse scattering maps per- 
formed along the lines perpendicular to the T vector (see 
text for the details). These lines are shown in solid on Fig. (2] 



structure factor. Summation is made over v possible dif- 
ferent orientations of defects (i^ > 1 if the symmetry of 
the defects is lower than that of the crystal). 

The equation ([T]) describes the deformations of an ideal 
anisotropic elastic medium with low concentration of de- 
fects and predicts an infinite increase of ^^(q) |g-j.o- The 



corresponding intensity of scattering is proportional to 
In real crystals one obviously could observe only 
finite intensities and one of the mechanisms that can sub- 
stantially smooth the profiles of Huang-type scattering 
near the ZC is lattice deformation screening. In the case 
of large concentration of defects the deformation fields 
produced by individual defects are spatially screened due 
to the infiuence defects in the neighborhood. This results 
in exponential cutoff factor for the deformation fields in 
the form 

u{r) - \e-'l''' 

where u{r) is the strain magnitude and Ts - the screen- 
ing length^. In the reciprocal space this fact is reflected 
by changing from l/g^ dependence to l/((7'^ + (ff). 

We account for this effect by modification of ((T|) in the 
form 

^»(q) = ^ NaCr^{ci)Vo.,ini (3) 

This way the scattering profies transform into 
Lorcntzians with parameter inverse proportional 
to the screening length r^. The experimental data show 
no sharp increase of diffuse intensity near the Bragg 
reflections (Fig. [71) and we attribute this fact to the 
presence of strong deformation screening in PMN. 

As it was mentioned before we consider the defects of 
cubic symmetry which are characterized by elastic defor- 
mation tensor of diagonal form Laji = SjiA. The defects 
of this kind expand or compress the surrounding lattice 
equivalently in all three main crystallographic directions, 
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the character (compression or expansion) is determined 
by the sign of constant A. In case of cubic defects Ih{Q) 
has a zero-intensity plane that is nearly normal to the 
vector T^^ and the intensity near the zone center can be 
approximated by expression const ■ cos{(j))/q where is 
the angle between T and q. This expression results in 
appearance of the thin waist near the zone center outside 
which the intensity is nearly zero (area A on Fig. Hljb)). 
Calculated 2-D distributions of DS intensity convoluted 
with experimental resolution function are presented in 
the bottom row of Fig. [5] One can clearly see that the 
three main peculiarities - the character of anisotropy, the 
presence of the waist and the satellite maxima are well 
reproduced by the model. 



IV. SATELLITE PEAKS 

In all studied Brillouin zones the diffuse scattering 
intensity contains maxima at finite q values. In high 
symmetry zones the peak positions lie on the line go- 
ing through the corresponding reciprocal lattice points. 
In zones of lower symmetry these positions are somehow 
shifted from that line. Equations ([3]) and ([2]) describe the 
intensity that monotonically decreases with increasing q 
value. In this section we show how finite experimental 
resolution can affect the scattering profiles resulting in 
the appearance of the observed peaks. In the case of 
elastic scattering the resolution function (RF) can be ap- 
proximated by 3-D gaussian^^: 



i?(q) = Ro 



o-ii{<]-<]o)f,/<y?,+{q-qo)±/<yj_ + {q-qo)v/<yv) 



(4) 



Here the values (q — (7o)|| ■•• {q^Qo)v a-re the components 
of vector (q — qo) along the vector T, along the in-plane 
and out-of-plane perpendiculars correspondingly. In our 
experiment the in-plane components of the RF were quite 
sharp as it can be deduced from the images of the Bragg 
reflections on Fig. [2][a) . In order to provide high neutron 
flux the vertical coUimations were highly relaxed which 
resulted in more broad resolution in out-of-plane (verti- 
cal) direction. This is illustrated on Fig. IH^b) by the 
i?(q) = -Ro/2 resolution ellipsoids stretched in (0 I) 
direction. 

The figuresHl^b) (2-D contour plots) and[5](l-D slices) 
represent the same distribution of the intensity calcu- 
lated using equations ([3]) and ([2]) in (3 0) Brilluin zone 
together with the resolution function. One can see from 
these figures that (a) the intensity decreases with increas- 
ing q value and (b) the distribution of intensity along 
directions is very narrow for small h values and is 
broadening on moving from the zone center along (h 
0). The first factor leads to the decrease of the detector 
count rate with increasing h value. The second factor 
on the other hand leads to its increase due to increasing 
of overlapping between the intensity distribution and the 
RF. Near the zone center the second factor is dominant 
which results in the increase of the count rate. This trend 
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FIG. 4: a) A set of modelled radial scans of diffuse scat- 
tering (without Bragg and background contributions) near 
(3 0) reflection with difi^erent resolution conditions. The 
thick line corresponds to the estimated resolution of current 
experiment, dashed line - to the estimated resolution of ex- 
periment from Ref. [l^. The uppermost line corresponds to 
the ideal resolution, the second upmost line corresponds to 
AQy = ^Qz ~ 0.011 A~^, for other lines the resolution 
conditions are denoted on the plot, b) The projections of 
resolution ellipsoids (thick lines) for the zone center (A), q—- 
0.125 (B) and g=0.125 (C) and isointensity contours of the 
Huang scattering with parameters described in text. Resolu- 
tion function and scattering intensity are normalized to unity, 
the thick lines correspond to the value 1=0.5. 



changes near h=Q.15 r.l.u. where the widths of the in- 
tensity distribution and the RF become comparable and 
the impact of the second factor diminishes. On further 
increase of h the detector count rate follows the decrease 
of the unconvoluted intensity (Fig. SJa)). The crossover 
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FIG. 5: The slices of 3-D model scattering intensity in (3 0) 
Brillouin zone determined by equations ((3]) and ([2]) along the 
vertical I axis with fixed value of fc=0 and different values of 
h described on the graph. Bold line represents the normalized 
resolution function. 
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FIG. 6: (Color online) Calculated diffuse scattering inten- 
sity in (110) Brillouin zone. Corresponding experimental data 
were reported in Refl2^. 
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between these two regimes is marked by satellite maxima 
at (7=0.15 r.l.u. 

The maxima do shift towards larger q values with in- 
creasing tj values. On Fig|4lja) the profiles for different 
cry values are presented including the conditions of our 
experiment (bold solid line) and experiment by Hiraka et 
al.-^S (dashed line). Their vertical resolution should be 
about 2.6 times better than our due to use of cold neu- 
trons with Ei — Ef = 4.5 meV. The resolution effect in 
this case is is still strong enough to produce the observed 
satellite maxima. The 2-D map calculated with these ex- 
perimental conditions for (1 1 0) BZ (Fig. [5]) is in very 
good qualitative agreement with the experiment. 

In all calculations presented in this paper we used the 
reciprocal screening length Qs — 0.2 which we found to 
be optimal for the description of the scattering profiles 
(FigEl). 



DISCUSSION 
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FIG. 7: The longitudinal (radial) scan of diffuse scattering in 
PMN at T=650 K in the vicinity of (3 0) Bragg reflection 
(circles), model fit (solid line) and the decomposition of fit 
into the Bragg, diffuse and background components (dashed 
lines). 



We should discuss the microscopic characteristics of 
PMN that can be responsible for the observed diffuse 
scattering. Relaxors are mixed crystals and are definitely 
expected to suffer from elastic deformations due to the 
intrinsic inhomogenity. Indeed the deformations mod- 
elled on the basis of known in advance elastic constants 
and the most simple cubic symmetry defects are shown 
to produce the scattering with the same shapes as ob- 
tained in the experiment. The role of defects most prob- 
ably play the inhomogenities of chemical composition. It 
is naturally to expect that the relative density of Mg'^+ 
and Nb^"*" ions is weakly fluctuating on the nanoscale. 
The areas characterized by excess of one type of these 



ions are expected to have different from average unit cell 
volume due to difference in effective ionic radii (0.72 Afor 
Mg^+ and 0.64 Afor Nb^+)^ and thus produce expan- 
sion or compression of the surrounding lattice. The areas 
with excess of one or another type of ions most proba- 
bly should be considered not as well-localized entities but 
as entities smoothly passing into each other. The cross- 
overs between regions with excesses of different ions take 
place on the nanoscale so it is naturally to assume that 
the concentration of defects is large enough to produce 
deformation screening needed to suppress the sharp in- 
crease of diffuse intensity near the zone centers. Possibly 
the chemically ordered regions also contribute to the lat- 
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tice deformations due to the different from the disordered 
matrix unit ceh volume. 

Also we address the problem of satellite maxima near 
the Bragg reflections. Typically the satellites are created 
when the new periodicity longer than the unit cell is cre- 
ated. The positions of the satellites are determined by 
the ratio of the "new" and "old" periods . In the case 
of high-temperature diffuse scattering in PMN the posi- 
tions of maxima vary between different BZs but remain 
on the nearly the same distance from BZs center. If we 
assume that the peaks originate from short range corre- 
lations of ionic displacements as it was suggested in Rcf. 
[2^ it seems to agree well with anomalously large width of 
the peaks but it also seems very difficult to construct a 
quantitative model predicting the positions of the peaks. 
The model proposed here together with the resolution ef- 
fect provide on the other hand a simple way of satellite 
peak description without any additional assumptions ex- 
cept the lattice deformation screening which by itself is 
quite expectable in relaxors. 

VI. CONCLUSION 

In this contribution we report the results of the de- 
tailed study of the high-temperature diffuse scattering 



in low-symmetry Brillouin zones in PMN and propose a 
model that accounts well for both the discovered scat- 
tering anisotropy and the presence of satellite maxima. 
The assumptions on which the model is based are simple 
and consistent with the data obtained earlier. To clarify 
the microscopic origin of the elastic lattice deformations 
additional measurements are necessary. 
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